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The confluence of quantum physics and biology is driving a new generation of quantum-based 
sensing and imaging technology capable of harnessing the power of quantum effects to provide tools 
to understand the fundamental processes of life. One of the most promising systems in this area is 
the nitrogen-vacancy centre in diamond - a natural spin qubit which remarkably has all the right 
attributes for nanoscale sensing in ambient biological conditions. Typically the nitrogen-vacancy 
qubits are fixed in tightly controlled/isolated experimental conditions. In this work quantum con- 
trol principles of nitrogen-vacancy magnetometry are developed for a randomly diffusing diamond 
nanocrystal. We find that the accumulation of geometric phases, due to the rotation of the nan- 
odiamond plays a crucial role in the application of a diffusing nanodiamond as a bio-label and 
magnetometer. Specifically, we show that a freely diffusing nanodiamond can offer real-time infor- 
mation about local magnetic fields and its own rotational behaviour, beyond continuous optically 
detected magnetic resonance monitoring, in parallel with operation as a fluorescent biomarker. 

PACS numbers: 03.65.Vf, 03.65.Yz, 42.50.Dv, 76.30.Mi, 05.40.Jc, 87.15.H- 



INTRODUCTION 

The diamond nitrogen-vacancy (NV) defect shows 
great promise as a biological imaging tool in two dis- 
tinct capacities. Firstly, its optical properties and bio- 
compatibility make it an excellent fluorescent biolabel 
[1-8]. It is extremely bright, highly photostable [9], and 
it can be imaged with spectacular spatial resolution [10] 
through sub-diffraction localization techniques like stim- 
ulated emission depletion [11]. Secondly, optical spin po- 
larisation and measurement and microwave control, the 
NV defect's electron spin can be used to detect magnetic 
fields with nT Hz -1 / 2 [12, 13] sensitivity, with potentially 
nano-metre spatial resolution [14, 15]. The prospect of 
using the NV centre for nanoscopic MRI has significant 
promise for a range of problems in biology, including: de- 
tecting operation of individual neuronal ion channels [16] 
and atom-scale magnetic resonance imaging [17, 18]. 

However, the more revolutionary application is the 
combination of the NV as a nano-scale biomarker and 
magnetometer. This would be a super-optical probe ca- 
pable of reporting motional information simultaneously 
with atomic level magnetic fields and their fluctuations 
over arbitrary long times. Such a system would be able 
to give the precise location of some attached ligand as 
it moves through a cell, as well as information about 
its molecular environment by detecting nanoscale mag- 
netic fields. However, nearly all demonstrations of NV 
magnetometry have assumed the NV spin quantization 
axis is effectively stationary with respect to the lab- frame 
control fields. In order to employ the NV centre as a 
combined nanoscale bio-mar ker/magnetometer one must 




Figure 1: (a) Energy level diagram of the NV center, (b) Ge- 
ometry. The magnetic field direction of the Rabi pulse defines 
the z-axis. z' is the instantaneous direction of the N-V axis, 
defined with respect to the lab-frame unprimed coordinate 
system by 8, (j>. 



fully understand the effect of random rotational motion 
on the nanomagnetometry and control protocols. In Ref. 
[19] NV nanodiamond nanomagnetometry in living cells 
was demonstrated. In these experiments, the measure- 
ments were conducted on a endocytosed nanodiamond 
which moderated the motional timescales, facilitating 
quantum measurement and allowing information about 
the orientation to be determined over long periods. In 
this paper we explore the quantum control of a rotation- 
ally diffusing NV system and find that the accumulation 
of geometric phases, due to rotation [20-25], is the criti- 
cal phenomenon in the application of the NV centre as a 
bio-marker/magnetometer. 
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THE DIAMOND NV SYSTEM 

The NV defect has a spin triplet ground state with a 
2.88 GHz zero-field splitting between the m — state 
and the degenerate m — ±1 states as shown in Fig. 1(a). 
Optical excitation at 532 nm can pump the defect into 
the m = state, and allows the population of the ground 
state to be read, since the m = state has a higher flu- 
orescence than the m — ±1 states. The effective Hamil- 
tonian of the NV ground state triplet in the presence of 
a magnetic field B [lab frame coordinate system defined 
by z, see Fig 1(b)], ignoring crystal asymmetries and hy- 
pcrfine effects, is 



H 



1 



-DSi 



jBS, 



(1) 



where 7 is the hyro- magnetic ratio of the NV (7 = g^s/H, 
g ss 2). The first term is the zero-field splitting of the NV 
system itself, where D is the zero-field splitting strength. 
It is this term which makes the crystal's orientation cru- 
cial. It defines a quantisation direction z' [see Fig. 1(b)] 
which lies along the axis connecting the nitrogen atom 
to its adjacent vacancy. The second term is the usual 
Zeeman splitting interaction with a magnetic field. 

Previously it has been shown that the NV centre can be 
used to measure DC and AC magnetic fields [12, 14] and 
fluctuating (FC) fields [26, 27] and rotation [19-21, 23- 
25]. These protocols work by optically polarizing the NV 
m = spin state. A 7r/2 microwave pulse then estab- 
lishes a coherent superposition of the m — and m = 1 
states, and the system evolves under the influence of any 
magnetic fields by the Zeeman effect. After an evolu- 
tion time, t, the rotation and strength of the magnetic 
field manifests itself as a phase, relative to m — 0, of the 
m = 1 state. In the adiabatic limit this phase is 



$ = $ G + 3> D = 



cos6»d0 + 7 / z'(t) ■ B(t)dt, 



(2) 



where $ G is the accumulated geometric phase, due to 
rotation along some path, P, defined by 9 and [see Fig. 
1(b)] [21, 23] and $ D is the accumulated dynamical phase 
due to the change in the projection of the magnetic field, 
B onto the NV axis, z! . A second ir/2 pulse converts 
this phase into a population difference, which is read out 
directly as a change in fluorescence. 

It is common to apply a static magnetic field to split 
the m = ±1 states, selecting only one of the m = ±1 
states by using a sufficiently weak microwave field res- 
onant with one of the transitions (Rabi field). Since 
a static field would cause different splitting depending 
on the crystal's orientation [19], we do not take this ap- 
proach. Instead, we assume that the degeneracy of the 
m = ±1 states is lifted by strain, which is usually the 
case for nanodiamonds [19]. 
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Figure 2: Important timescales as a function of crystal ra- 
dius, assuming the crystal is immersed in room temperature 
water (T = 300 K) with a viscosity of 77 = 10~ 3 Pa s. Red 
dashed curve: characteristic rotational diffusion timescale, 
1/fcd- Blue dotted curve: population mixing timescale, t m , 
for non adiabatic motion. Green dashed-dotted curve: angu- 
lar velocity damping timescale, td- Also shown are the typical 
timescales associated with the NV spin state: homogeneous 
broadening (T2), inhomogeneous broadening (TJ) and zero 
field splitting (1/D). 

MODEL OF DIFFUSION 



Motion of the crystal 

A spherical crystal is considered, whose rotation fol- 
lows Brownian motion, according to the Langevin equa- 
tion 



Ibj = -7 d w + r(t), 



(3) 



where / is the crystal's moment of intertia, u> is the crys- 
tal's angular velocity, r is a stochastic delta function 
correlated torque, and 7<j is a drag coefficient given by 
Stokes' law, 7^ = 87r?7r 3 , where r is the crystal's radius 
and r\ is the fluid's viscosity. It can be imagined that 
the crystal rotates in a 'random walk' of discrete small 
angle steps, whose directions are uncorrelated with one 
another. This rotational motion is characterised by a sin- 
gle diffusion coefficient kd which, from Eq. (3) and the 
equipartition theorem, is given by 



kd = 



k B T 



(4) 



where k B T is the thermal energy. Thus, in room- 
temperature water, r = 5 nm [r = 50 nm] gives kd ~ 1 
rad 2 //xs [kd ~ 1 rad 2 /ms], see red dashed curve in Fig. 2. 



Ensemble averaged measurements 

In practice, and even in theory, a single measurement, 
of a quantum system, contains very little information. 
The quantity which is measured is a time averaged signal 
from a large ensemble of individual measurement runs. 
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Environmental perturbations cause subtly different quan- 
tum evolution in each member of the ensemble, leading 
to a loss of information in an ensemble-averaged signal 
after a certain evolution time. In the case of the rota- 
tionally diffusing nanodiamond, the ensemble's members 
differ more dramatically. In addition to its own micro- 
scopic environment each member of the ensemble has its 
own starting orientation and rotational trajectory. 

A series of sequential measurements is made which may 
involve several NV centres in a single crystal. If the time 
over which measurements are made is sufficiently small 
compared with the rotational diffusion time, 1/kd, then 
the situation would be comparable to the case of a static 
crystal. In this work the opposite limit is considered and 
it is assumed that sufficiently many measurements are 
made that the final signal is an average over the full the- 
oretical ensemble: all initial orientations and rotational 
trajectories are accounted for. Experimentally this aver- 
age can be constructed via the repeated measurement of 
the same NV centre or over a collection NV centres [19]. 

For such a scenario consider a continuously applied mi- 
crowave field, of amplitude Br, tuned to the ijr, — > 
^(m-i) ^ rans j^j on _ This will cause Rabi oscillations with 
a frequency which depends on the angle, 0, between the 
NV axis z' , and z, the oscillation direction of the mi- 
crowave magnetic field. Initialised into ip^\ the popu- 
lation of ipf) for a single NV as a function of time is 
then 

P = [l + cos(n fl tsin0)]/2, (5) 

where £Ir = g^BBuJh is the Rabi frequency when 
6 = 7r/2. Assuming £Ir 3> 1/kd the ensemble-averaged 
signal is an average over all the possible different Rabi 
frequencies: 

i r 

S (t) = - / d9 sin 6 [1 + cos (n R t sin 6)} 

= J [2 + 7rff_i (n R t)] , (6) 

where H a is the Struve function. So(t) (solid curve) is 
shown in Fig. 3(a), along with the equivalent curve for a 
static crystal (dashed curve), with 6 = ir/2. For an en- 
semble averaged measurement the signal reaches its first 
minimum at tfl R w l.lfwr at which point So ~ 0.14. 

This result enables the determination of the optimal 
microwave pulse time to achieve a 7r/2 rotation. The re- 
quired time is simply read off from the curve. In our 
system, there is a certain amount of freedom in decid- 
ing how long a 7r/2 pulse should be. We could take it 
to be ijr/2 = 7r/(2£2fl), which would produce a n/2 rota- 
tion for an optimally aligned (8 = 7r/2) NV centre, and 
a smaller rotation otherwise. Alternatively, we could use 
= l-167r/ (251/;.) which, based on Fig. 3(a), would 
optimise the initial amplitude of a Ramsey-type experi- 
ment. One of the strengths of Ramsey- type experiments 




Figure 3: (a) Rabi signal averaged over the ensemble of pos- 
sible crystal orientations (solid black curve) and for optimal 
crystal orientation (9 = n/2) (dashed green curve), (b) Signal 
envelope showing exponential decoherence due to geometric 
phase accumulation, for a = 1 (solid curve) and a = 1.16 
(dashed green curve). 

is their robustness to the ir/2 pulse time chosen. We 
do not, therefore, expect this decision to be an impor- 
tant one (as simulations will confirm). For now we keep 
our discussion general, introducing a parameter, a, to de- 
scribe the pulse time. As such when we refer to a '7r/2' 
pulse we mean a pulse of length f_/ 2 = omj ) , where 
a « 1, which produces a spin rotation of = a(7r/2)sin#. 

Now consider a Ramsey- type pulse sequence. The NV 
centre, starting with z' at some polar angle d\ to the z- 
axis, is first optically pumped into %jy^} . A l 7r/2' pulse 

is then applied, tuned to the tp^P — > tpj, transition, 
producing a spin rotation of ©i = a(7r/2) sin#i. The 
system is allowed to evolve for time r, during which 
time the crystal rotates through some trajectory and the 
tp^/ state develops a phase $, due to rotation and the 
presence of magnetic fields. A final ! 7r/2' pulse is ap- 
plied. This final pulse may be in phase with the initial 
pulse, or out of phase by tt. It produces a rotation of 
±0 2 = ±a(7r/2) sin 2 where 02 is the final angle between 
z' and z and the + (— ) refers to the pulse being in (out) 
of phase with the original pulse. The final population Pq 
of the ij/i) state can then be written down explicitly as 

P = ^ [1 + cos 9i cos 2 ± sin 0i sin 2 cos $] . (7) 

The final, normalised, signal is the ensemble average of 
this value S(t) = (Po), where the ensemble average is 
taken over all possible starting orientations and rota- 
tional trajectories. The upper and lower bounds of the 
signal envelope are produced by the two curves implied by 
the ± term of Eq. (7), provided the mean phase vanishes, 
(<&) = (true for the cases considered in this work). This 
work primarily focuses on the behavior (typically decay) 
of the signal envelope, due to rotation and the presence of 
magnetic fields, rather than the oscillations themselves. 
Physically, the two curves correspond to the phase (0 or 
7r respectively) of the final l ^^/2 , pulse. 
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DEPHASING DUE TO GEOMETRIC PHASE 
ACCUMULATION 

To see the influence of the geometric phase, we con- 
sider how the spin state of the rotating crystal evolves 
in the absence of external fields. As the crystal rotates, 
the direction of z' changes, so that Eq. (1) becomes a 
time-varying Hamiltonian. With respect to the eigen- 
sates of spin projection along the fixed z axis, the zero- 
field Hamiltonian is: 



- -D") 2 



Dh 



V2 



l_ cos t> sin t 

sin 2 6 



' cos 9 sin 8 
V2 



cos 8 sin 6 „„„2 a i sin 2 
■s/2 + 2 



(8) 



where 9 and 4> are time dependent. If the rate of rotation 
is small compared with 2.88 GHz the evolution will be 
adiabatic, producing no changes in populations of the 
spin sublevels (with respect to the z' quantisation axis). 
Each sublevel, however, develops a geometric phase 



$m = 







10 


dt 



(9) 



There is a gauge degree of freedom in defining the phase 
of the eigenstates [21, 23]. We choose to define the eigen- 
states with respect to the fixed z basis as 

1 „~i4> a\ 



cos 2 \ 

2^ sin 2 | 
sin 2 | 

cos 2 | 



(0) 



V2 



COS 

e** sin( 



This convenient choice of gauge means that an explicit 
phase factor between sublevels corresponds to the abso- 
lute phase of a microwave field linearly polarized along 
the z axis. The relative geometric phase accumulation, 
between if& and iffi , due to the crystal's rotation then 
evolves according to $ = (f> cos 9, where the phase evolu- 
tion of each state is given by $ m = m<£> [20, 21, 23-25]. 

Given a particular initial orientation 9\ the probability 
density of a particular final orientation 9 2 and geomet- 
ric phase $ after time t is p{9 2 , $; t\8\). The ensemble- 
averaged signal can then be expressed as the integral of 
the final ground state population over all possible 9\, 62 
and weighted by the probability distribution function 
and the probability of the NV starting with orientation 
9x: 



S{t) = - J d8 1 d9 2 d^>sm9 1 p(9 2 ,^;t\9 1 )P > 



92,$), 

(10) 




Figure 4: Rotation rate and DC magnetic field sensitivities, 
for an assumed fluorescence collection efficiency of a = 0.01 
and a density of NV centres in the crystal of 10 16 cm -3 . Solid 
black curve, right axis: DC magnetic field sensitivity. Red 
dashed curve, left axis: rotation rate sensitivity from geomet- 
ric phase decoherence. Blue dotted curve, left axis: rotation 
rate sensitivity from non-adiabatic population mixing. 

where P o (0i,0 2 ,<fr) is given by Eq. (7). The Fokker- 
Planck equation for the evolution of p (9 2 , t\9\) is 



~p(0 2 ,<M0i) 

kd dt 



d 

W 2 



sin 9 2 



d 
W 2 



P (9 2 ,^t\8 1 ) 



sin t) 2 



cot 2 9 2 —p{9 2 ,<5>;t\9 1 ) 



(11) 



The first term describes the probabilistic rotation of the 
crystal, and the second term describes phase evolution, 
which dctcrministically depends on the probabilistic evo- 
lution of azimuthal angle (j). The initial condition is 



p(9 2 ,$,t = 0\9i) = 8{$)5(e 2 -9 1 ) ) 



(12) 



and the boundary conditions are periodic in 9 and $. 

Numerically solving Eq. (11) gives the solution shown 
in Fig. 3(b), demonstrating that the decaying envelope is 
relatively insensitive to the l 7r/2' pulse duration, a. There 
is approximately exponential decay of the signal envelope 
over the timescale of the crystal's rotation, with decay 
time « 0.90/^. The diffusion rate, kd, thus sets an up- 
per bound on the coherence time for the crystal. In the 
absence of any external fields, measurement of this deco- 
herence rate allows measurement of kd with shot noise- 
limited sensitivity AkdVT — where is the 

overall decoherence time (l/ra ~ ^/T 2 + kd), a ps 0.01 is 
a parameter which accounts for the non-zero fluorescence 
of the m = ±1 states, and finite collection efficiency [17]. 
N is the number of NV centres in the crystal, along all 
four crystallographic axes, and T is the total measure- 
ment time. Figure 4 (red dashed curve) shows the sen- 
sitivity as a function of crystal radius, with an assumed 
NV density of 10 16 cm- 3 . 



DEPHASING DUE TO POPULATION MIXING 
WITH NON-ADIABATIC EVOLUTION 

Up to this point we have been working within the adi- 
abatic approximation. In the context of the NV centre, 



5 



the adiabatic criterion is 



6 2 sin 2 i 



< D 1 



(13) 



The angular velocity of the nanodiamond is a conse- 
quence of thermal motion. Using the equipartition theo- 
rem, the adiabatic criterion for the rotationally diffusing 
nanodiamond is 



I 



= 15fcd, 



(14) 



8jrpr 3 

where p is the density of diamond. At room temperature 
this criterion will be satisfied if the crystal radius is much 
greater than 0.2 nm [see Fig. 2, dashed red curve]. 

An additional adiabatic criterion, often overlooked, is 
that the state of the NV centre changes slowly compared 
to 1/D. The timescale on which the NV state changes 
is the angular velocity damping time, td = I /id — 
pr 2 /(15r/). In Fig. 2(a) t d (green dashed-dotted curve) 
is plotted as a function of the radius of the crystal, for a 
crystal in room temperature water. From this it can be 
seen that td is smaller than 1/D when the crystal radius 
is smaller than 20 nm. 

The geometric phase accumulation comes from a ze- 
roth order perturbative expansion of the time evolution 
operator in the ratio of the rotation rate to the 2.88 
GHz zero field splitting. Taking the expansion to higher 
orders enables the population mixing that comes from 
nearly adiabatic evolution to be evaluated. To evaluate 
the properties of an ensemble measurement consider an 
arbitrary initial state of the NV centre, where each state 
has a probability weighting of P m and phase <I> m . As 
shown in Appendix A the ensemble average population 
mixing rate is 

<*> = rr'SW (15) 



l 



^ 2 • 



The above provides an expression for the population mix- 
ing rate in terms of physical parameters: the zero-field 
splitting frequency D, the rotational diffusion rate kd and 
the angular velocity damping time td- 

To obtain the above results several assumptions have 
been made. Most significantly, the diffusive rotation of 
the nanodiamond has been modeled as a series of dis- 
crete steps, using physical arguments to infer the appro- 
priate statistical behaviour of the steps' sizes and dura- 
tions. To check the validity of these approximations a nu- 
merical Monte Carlo simulation has also been employed. 
The simulation models the Langevin equation (3), us- 
ing a torque which fluctuates on a finite but very small 
timescale, to produce a set of rotational trajectories. For 
each trajectory the Schrodinger equation is solved numer- 
ically. The results, for an NV centre which starts in 
are plotted in Fig. 5(a). They show excellent agreement 
with our analytical expression: 

S m (t) = - + ^exp(-t/t m ). (16) 
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Figure 5: (a) Population mixing due to non-adiabatic evolu- 
tion: expected fluorescence signals as a function of free evo- 
lution time, relative to rotational difusion time, given initial 
occupation of jO) . The dashed curves follow the theoretical ex- 
ponential decay of Eq. (16) while the solid curves are based 
on a Monte Carlo simulation of the Langevin equation fol- 
lowed by numerical solution of Schrodinger's equation. Blue 
tdD = 0.1, red: tdD = 1, green: tdD — 2. (b) Decoherence 
curves for several DC field strengths and orientations. Solid: 
no field. Dashed: B = Wkd/fB. Dotted: B = 10k d /yx. (c) 
Decoherence rate relative to kd for DC magnetic fields in z 
(solid) and x (dashed) directions, (d) Decoherence signal as a 
function of field strength for a spin-echo pulse sequence with 
an AC field in the limit of slow rotation compared to free 
evolution time. 

where t rn = l/(3(fc)) is the population mixing time. 

The implication of this non-adiabatic behaviour is 
more than merely a correction to the behaviour ex- 
pected from geometric phase accumulation. Measuring 
the population mixing rate (effectively measuring Ti) 
provides another way to gauge the rotational diffusion 
rate. More over the measurement can be achieved with- 
out microwave fields. It is enough to simply initialise the 
NV into the ip^ state and observe the decay of the flu- 
orescence rate as a function of free evolution time. The 
blue dotted curves in Figs. 2 and 4 show the population 
mixing time, and associated sensitivity, respectively, as a 
function of the crystal radii. 



STATIC AND FLUCTUATING MAGNETIC 
FIELDS 

We now investigate the properties of rotationally dif- 
fusing NV centres in the presence of static (DC) , fluctuat- 
ing (FC) and oscillating (AC) magnetic fields. Magnetic 
fields enter our model by advancing the phase of the NV 
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centre's spin sublevels. Under the secular approximation 
Eq. (1) becomes 

^DS% +-f(B z cos0 + B x sm8cos(f)) 7 (17) 

it is only the field component parallel to z' (the NV axis) 
which affects the phase. 



field, the magnetic field sensitivity is 
ABVT; 



(19) 



where 1/r^ s=s T 2 + kd- The solid black curve in Fig. 4 
shows the DC magnetic field sensitivity as a function of 
the crystal radii in room temperature water. 



Static magnetic fields 



Fluctuating magnetic fields 



Since the orientation of z' changes randomly with time, 
a static field effectively behaves as a field which fluctu- 
ates on the timescale of rotational diffusion 1/kd [27]. 
To demonstrate this consider a static magnetic field with 
components B z and B x along the z and x axes (without 
loss of generality, B y = 0). Unlike the geometric phase, 
the phase evolution due to a magnetic field depends ex- 
plicitly on the azimuthal angle <j> as well as the polar 
angle 6 of z'. The Fokker-Planck equation is then: 



sin 02 ggr 



sin 62 



- cos 2 + ^ sin 9 2 cos 4>) ^P{02,4>, *, t|*i)- 

(18) 

The resulting signal envelope [Eq. (10)] is plotted in Fig. 
5(b) and the decoherence time (defined here as the time 
taken for the signal envelope to reach 1/e of its initial 
value) is plotted in Fig. 5(c). 

When jB 3> kd, the shape of the dephasing signal is 
close to Gaussian, the characteristic shape for slowly fluc- 
tuating fields. In this regime the static field effectively 
fluctuates on the timescales of 1 /kd due to the rotation of 
the crystal. The crossover between Gaussian and expo- 
nential decay occurs on the timescale of 1 /kd, but the ex- 
ponential behaviour is hidden since the geometric phase 
accumulation causes dephasing to occur on at least this 
timescale. The decay is significantly stronger for fields 
in the x (or y) directions than in the z direction. This 
is because the largest Rabi rotations are experienced by 
the ensemble members for whom 8 w ir/2 which are most 
sensitive to fields in the x-y plane. 

For weak fields, the sensitivity approaches zero, which 
can be seen in the vanishing gradient for small B in Fig. 
5(c). To sense small DC fields, a strong static magnetic 
field could be applied to bring the field sensitivity to 
within an appropriate range. The nonlinearity of field 
sensitivity also allows the possibility of vector magnetom- 
etry. By applying a DC field in one direction, only DC 
fields in that direction would be detected, allowing both 
the magnitude and the direction of DC magnetic fields 
to be imaged. Given the optimal strength of applied DC 



Many biological systems display behaviour that is nei- 
ther static nor periodic, but stochastically fluctuating. 
Although the details may vary, a fluctuating field can 
be characterised by a correlation time t c , and a mean 
squared field strength (B 2 ) [27]. For the tumbling NV 
system, a magnetic field correlation time larger than the 
rotational diffusion time will mean that the NV experi- 
ences a field that is essentially static over its coherence 
time, although statistically distributed over many runs. 
In this case, the effect will be very similar to that for a 
DC field. 

For the case of a fluctuating field with correlation time 
t c <C 1/kd and mean squared field strength {B 2 ) = 
B 2 FC <C l/(7i c ) 2 which is isotropic, the Fokker-Planck 
equation is 



1 dp(0 2 ,M6»i) 



dt 



d_ 
W 2 



sin 2 



8 f P (0 2 ,<P;t\0 1 ) 



d0 2 



sin 2 



+ cot 2 2 + 



i 2 B 2 FC t c \ d 2 p(# 2 ,<M|0i) 

6k d 



<9$2 



(20) 



which is similar to that for the geometric phase alone 
(11). As such the fluctuating field produces exponential 
decay of the signal, on a timescale of l/{^B FC y/T c ) 2 . The 
implication of this result is that the decoherence rate is 
linear in Bp C t c . 



Oscillating magnetic fields 

One application of a nanoscale magnetometer is to de- 
tect magnetic resonance from a small number of spins, 
which can be achieved by measuring the AC field pro- 
duced by their effective precession, or by driving them 
via an external microwave field. Applying a spin-echo 
pulse sequence — ^- — 7r — — effectively rectifies 
an AC magnetic field of period 2r, allowing its ampli- 
tude to be detected while masking fields oscillating at 
different frequencies. A spin-echo pulse sequence has the 
additional advantage of dramatically increasing the de- 
phasing time of an NV magnetometer, by refocusing the 
divergent phases accrued due to slowly fluctuating mag- 
netic fields. 

In detecting magnetic resonance, a strong DC field is 
usually applied. A tumbling crystal will see such a field 
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as if it were fluctuating on the timescale of the crystal's 
rotation, causing dephasing on a timescale, T n7r , given by: 

T - ~ fej ' (21) 

where n is the number of 7r-pulses employed in a conca- 
tentated spin-echo pulse sequence. 

Taking the limit in which very little rotation occurs 
on the timescale of free evolution. The problem is then 
a matter of simply averaging over the possible crystal 
orientations and initial phases of the AC field, rather 
than solving a Fokker-Planck equation. Consider a spin- 
echo pulse sequence during which the crystal is oriented 
at some 9 and the free evolution time r matches the AC 
fields period. The final population of ip^P is given by 

P (0, Afield) = ^ + - cos 2 9i cos 2 9 2 

± ^sin 2 ei(l-cos 2 e 2 )cos(2$ AC ) 

(22) 

where 0i = a(n/2) sin# and 62 = a7rsin6> are the angles 
through which the spin is rotated during the tt/2 pulses 
and 7r-pulse respectively. For an AC field oscillating along 
the z-axis, the phase acquired by the NV centre before 
the 7r pulse is <&ac = ( t /2)7-Bac cos0cos<j>fi c id, where 
Afield is the phase of the AC field at the start of the 
period, and the phase acquired after the n pulse is — $ao 
In the limit that the strength of the AC field to be de- 
tected is much smaller than the applied DC field strength, 
it can be assumed that the crystal is stationary over the 
timescale of free evolution, although an average over the 
full ensemble of possible orientations still must be per- 
formed. To find the ensemble-averaged signal, we in- 
tegrate over all phases of the AC field, <&fi e id, and all 
orientations, 0, weighted by | sin 9: 

^(^f^^. (23) 

Figure 5(d) shows the signal envelope due to an AC field. 
As with DC fields, the optimum AC field strength is such 
that the dephasing due to the AC field matches the de- 
phasing rate from other sources. The optimal sensitivity 
to an AC field is the same as that given in (19), ex- 
cept with l/Td ~ 1/T 2 + l/r nv + kd- Using a concate- 
nated pulse sequence can increase the decoherence time. 
Concatenation also reduces the bandwidth of AC field 
detection, which is important if different AC frequency 
components are to be distinguished. 

CONCLUSION 



sensitivity comparable to that of a fixed crystal magne- 
tometer. We have also shown that, through the accu- 
mulation of a geometric phase, the fluorescence signal 
contains information about the crystal's rotational dif- 
fusion rate. This information can be gathered without 
applying any microwave control sequences, and is thus 
achievable with little modification to the existing experi- 
mental setups that use fluorescent colloidal nanocrystals 
as biological markers. 

Before the tumbling nanocrystal system can be used 
in biological applications, further theoretical and exper- 
imental work will be required. It will help to under- 
stand anisotropic rotation, and to investigate the effect of 
the fields produced by a specific biological process rather 
than the idealised fields considered here. The most chal- 
lenging experimental task will be finding protocols to ei- 
ther track or control rapid spatial movement of the crys- 
tal, while still achieving satisfactory microwave control 
of the quantum system within [19]. 
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APPENDIX A 

Ensemble average population mixing rate due to 
non-adiabatic rotation 

To evaluate the properties of an ensemble measurement 
consider an arbitrary initial state of the NV centre, where 
each state has a probability weighting of P m and phase 
$ m . The populations of the spin sublcvels after some 
time interval St is: 

P[ = Pi- kStP 1 + k5tP Q 

+ 2^/P 1 P VkSicos(<5> - $1 +H) (24) 
Pa = Pa- 2kStP a + kStPi + kStP^ 

- 2v / PiP v / A^tcos($o - $1 +H) 

- 2 v /P_iP \ / ^cos($ -$-i +£</>) (25) 
P'_ x = P_! - kStP^ + k6tP 

+ 2 v /P_ 1 P \/Mtcos ($ - $_! - 84>) , (26) 

where 



We have shown that a freely diffusing diamond 
nanocrystal can act as a sensitive magnetometer with 



(27) 
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Modeling the crystal's Brownian rotation as a random 
walk of small uncorrelated angular steps of duration St 
and angular change 50 on timescales much longer than 
St the ensemble average of P a can be determined: 



d(Pi) 
dt 

d(P ) 
dt 
d(P-l) 
dt 



-(fc)(P 1 ) + (fc)(P_ 1 ) (28) 
(k)(P 1 ) + (k)(P- 1 )-2(k)(P ) (29) 
-(fc)(P_ 1 ) + (fc)(P 1 ), (30) 



To find the ensemble average of k we need to know how 
St is distributed. If we multiply the Langevin equation 
(3) by u>(0) and take the ensemble average we find that 
the correlation function of the angular velocity decays 
exponentially over time, 



(u(t)-u(0)) = (u(0)-u,(0))e- t ^, 



(31) 



where we have defined the damping time, t d = I/"fd- In 
order to retain this correlation function in our discrete 
random walk approximation, the duration of a given step, 
St, must be exponentially distributed. At any given in- 
stant, the fraction of the ensemble that is experiencing 
a step of duration St is then given by the probability 
density: 



p(St) = 



-st/t d 



(32) 



We do not need to know the exact distribution of angu- 
lar velocities, merely that ((<56>) 2 /(^) 2 ) = 2/3(u> • u>) = 
2kd/td- The ensemble average population mixing rate is 
then given by 



(k) = 



J_/S(P 
D 2 \ St 2 

2k d 
l+t 2 D 2 ' 



(1 - cos(StD)) 



(33) 
(34) 
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